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Abstract 

The motion of a vortex filament in superfluid '^He is considered by using the Hall-Vinen 
([S], [ID]) phenomenological model for the scattering process between the vortex and thermal 
excitations in liquid ^He. The Hall-Vinen equations are formulated first in the intrinsic 
geometric parameter space to obtain insights into the physical implications of one aspect 
of the friction term, associated with the friction coefficient a, in the Hall-Vinen notation, 
as well as a proper rationale for the previous neglect of the other aspect of the friction 
term, associated with the friction coefficient a' . This development also serves to highlight 
the difficulties arising in making further progress on this route. A reformulation of the Hall- 
Vinen equation in the extrinsic vortex filament coordinate space is then given. This is shown 
to provide a useful alternative approach in this regard. Though the friction term associated 
with a, for very small a, has little capacity to make significant contribution to the vortex 
motion in a quantitative way, it is shown to be able to change the vortex motion aspects like 
the vortex kink characteristics in a qualitative way. This becomes possible because of the 
ability of the friction term to play the dual roles of driving force and drag force. 
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1. Introduction 



The Landau [1] model of superfluid ^He considers the superfluid below the lambda point as 
an inviscid, irrotational "background fluid" with thermal excitations, the phonons and rotons 
moving upon that background. These excitations constitute the normal fluid which would 
not interact with the superfluid except for the presence of vortice^. The vortices scatter 
the thermal excitations when there is a relative velocity between them, thus generating the 
so-called "mutual friction" (Feynman ^). The existence of mutual frictionH was confirmed 
by the experiment on the attenuation of second sound in uniformly rotating liquid ^He (Hall 
and Vinen [H], [ID])- These experiments, if the plausible assumption is made that a thermal 
excitation can exchange momentum only in a direction perpendicular to the scattering vortex 
line, then indicated (Vinen [11]) that the friction force will have a constant and large value 
for all directions of the relative velocity perpendicular to the rotation axis (with which the 
vortex lines are aligned) but will be zero for directions parallel to the latter. On the other 
hand, vortices in liquid ^He have a core radius of the order of the quantum coherence length 
(1° Angstrom), so the detailed core physics is not very relevant to the dynamical effects 
of the vortex away from the core. As a result, barring vortex reconnection events (which 
involve sharp distortions of vortex lines (Paoletti et al. [12]) and the concomitant generation 
of Kelvin waves associated with helical displacements of the vortex cores (Yepez et al. [13])), 
vortices in superfluid ^He behave like classical vortex filaments, the only difference being that 
their circulations and core radii exhibit quantum mechanical features. Indeed, the pioneering 
numerical simulations of Schwarz [Hj and [ihi, which provided considerable insight into 
superfluid vortex dynamics, are based on the idea that, except on very short length scales, 
quantized vortices can be regarded as vortex filaments moving according to classical fluid 
dynamics, with the inclusion of the mutual friction force. 

The leading-order behavior of the vortex-induced flow velocity formula in classical fluid 
dynamics is given by the so-called local induction approximation (LIA) (Da Rios [TB] , Arms 
and Hama [T7] ) which resolves the singularity due to the neglect of the finite vortex core size 
by an asymptotic calculation. Using the LIA, Da Rios [16] and Betchov [18] derived a set of 
coupled equations governing the inextensional motion of a vortex filament in an irrotational 
fluid in terms of time evolution of its intrinsic geometric parameters - curvature and torsion. 
Hasimoto [19] showed that Da Rios-Betchov equations can be elegantly combined to give 
a nonlinear Schrodinger equation. The single-soliton solution of this equation (Zakharov 
and Shabat [20]) provides a description of an isolated loop of helical twisting motion along 
the vortex line. The LIA is, however, hampered by the fact that the motion of the vortex 
filament is assumed to be governed solely by the local features on the filament, so distant 



^Vorticity in superfluid *He is confined to vortices which, as Onsager [2] suggested, are linear topological defects with the 
superfluid density vanishing at the vortex core. This allows liquid *He to have the angular momentum as in a solid body rotation 
situation while maintaining irrotationality, as per Landau's [l] theory, over almost the whole volume. This, on the other hand, 
as Onsager [2] pointed out, leads to the result that the circulation around a vortex line is quantized, which was confirmed in 
an ingenious experiment due to Vinen [3] . Rayfield and Reif [4] also gave direct evidence for the existence of quantized vortices 
in superfiuid *He via charge carrying vortex rings having one quantum of circulation. Direct observation of quantized vortices 
has been accomplished in superfuidity manifested in dilute Bose-Einstein condensates produced by laser cooling in laboratory 
experiments (Pethick and Smith 5 ). Direct observation of vortex cores has also been accomplished recently (Bewley et al. [g]) 
by using small solid hydrogen particles as traces in liquid *He. 

■^Our understanding of the microscopic origin of the physical processes underlying the mutual friction is not good (Donnelly 
[5]) so an adequate theory of the fundamental roton-vortex scattering process in superfiuid ''He is not at hand yet. 



2 



parts of tfie filament need to remain sufficiently separated during the motion. The large- 
amplitude solutions in LIA suffer from the violation of this premise as in a self-interaction 
of the vortex filament. 

Interestingly, rotating liquid "^He appears to be a better system than ordinary fluids for 
the application of LIA, because, 

• the cores of vortices in superfiuid "^He are only few angstroms in diameter (hence validat- 
ing the asymptotic evaluation implicit in the LIA which holds in the limit of vanishingly 
small vortex core size); 

• interactions between different segments of a vortex filament are important only when 
the distance between them is of the order of few angstroms ([S], [T5]). 

Upon including the effect of the frictional force exerted by the normal fluid on a vortex 
line, the self-induced velocity of the vortex line in the reference frame moving with the 
superfiuid according to the LIA is given by the Hall-Vinen equation (|9], [TO] ) 

V = 'jni X fi + at X (U — "yni x n) — a'i x [t x (U — 'jni x fi)] . (1) 

Here, U is the normal fluid velocity taken to be constant in space and time and prescribecH 
([14]. [T5]). K is the average curvature, and t and n are unit tangent and unit normal 
vectors respectively, to the vortex filament, and 7 = rin{c/ nao), where F is the quantum 
of circulation, c is a constant of order 1 and Oq ~ 1 • 3 x 10~^ cm is the effective core 
radius of the filament, a and a' are the friction coefficients which are usually found to be 
small (except near the A-point) so the short-term vortex motion appears to be only weakly 
affected by the friction. However, the friction term associated with a (which reflects the fact 
that a thermal excitation can exchange momentum only in a direction perpendicular to the 
scattering vortex line) plays the dual roles of driving force and drag force ([H], [IS]). It can 
therefore lead to both growth and decay of the vortex line length and hence can produce 
important qualitative effects. On the other hand, the friction term associated with a' arises 
partly from the asymmetry in the fundamental roton- vortex scattering and partly from the 
Magnus effect and is usually dropped under the pretext that a > au (Vinen and Nimela 
[21]). Nevertheless, determination of the vortex motion from the Hall-Vinen equation (1) is 
a highly non-trivial task and numerical simulations ([H], [l5]) have been so far essentially 
the only method of investigation. In this paper, we will first do Hasimoto [19] formulation of 
this problem in the intrinsic geometric parameter space to obtain insights into the underlying 
physical features as well as highlight difficulties involved in the Hasimoto [19] formulation 
route. We will then show that a reformulation of the Hall-Vinen equation (1) in the extrinsic 
vortex filament coordinate space provides a useful alternative approach in this regard - it 
provides insight into the fundamental importance of the friction term associated with a as 
well as a proper rationale for the previous neglect of the friction term associated with a' 



^So this formulation is kinomatical in nature - the effect of the vortices on the normal fluid is ignored. Notwithstanding 
this limitation, this formulation turns out to be useful in shedding light on the physical mechanisms involved. In actuality, U 
should be determined as part of the solution, once the back reaction of the vortices on the normal fluid has been accounted for. 

^Actual numerical values of a and a' at a couple of typical temperatures taken from table 1 of Schwarz |14| are - 

T{K) = 1.0 : a = 0.005, a' = 0.003 
T{K) = 1.5 : a = 0.073, a' = 0.018. 
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(as in the numerical simulations ([H], [15])). The friction term associated with a, even for 
very small a, is shown to be able to change the vortex motion aspects like the vortex kink 
characteristics in a qualitative way. 



2. Formulation of the Hall-Vinen Equation in the Intrinsic Geo- 
metric Parameter Space 

Let us first note that equation (1) can be rewritten as 

V = (1 — a') 7Kb + at X U + a'^nn (2) 

where b is the unit binormal vector to the vortex filament. Observe that the friction term 
associated with a' can be eliminated by renormalizing the vortex strength 7. 
We take the normal fluid velocity to be 

U = t/sfi + U^h (3) 

and allow the normal fluid velocity to be a function of the arc length s. Equation (2) then 
becomes 

V = 7Kb + a ( — t/sfi + f/2b j + a7Kn. (4) 

Note the Frenet-Serret formulae for the underlying differential geometry, 

x' = t, t' = Kn, n' = rb — Kt, b' = — rn (5) 

where r is the torsion and primes denote differentiation with respect to the arc length s. (5) 
leads to ^ 

fi + zb I = —ir I fi + zb I — Kt 



which suggests we introduce ([18]) 

N = (^fi + ib) e*^"(')'^^ i) = K(s)e*^"(')^^ (6) 

Note the following relations, 

N ■ t = 0, N ■ N = 0, N ■ N = 2 (7) 

where the bar overhead denotes the complex conjugate of the quantity in question. 
(6) leads to 

N' = -ijji (8) 

and 

t' = Re {ipN) = i (^N + V'N) . (9) 
On the other hand, we have from equations (4), (5) and (6), 

t = 7 ^K'b — Krnj + 07 ^Krb + k'tlJ — ar ^f/2fi + U^hj + aKU^i — a'^K^t 
= -fRe (#'N) + a J Re (tp'N) - a-f\ij\H + aRe (itpV) i + aRe {iV'N) 



i- {iIj'N - V^'N) + — (^'N + tp'N) - ajltplH + i- {tpV - ipV) i + i- {V'N - V'N) 
^ ^ 2 2 ^^^^ 
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where the dot overhead denotes differentiation with respect to time t, and 

V = {U2 + iUs)e'^^^'^'^'. (11) 

Let 

N = aN + (12) 

and on using (7) and (10), 

a + a = - (n ■ N + N ■ n] = -— (N ■ N) = or a = iR (13) 

and 

/3 = N . t = -N ■ t = - (i7 + a7) V-' - iaV (14) 

where R is some real-valued function. 

Substituting (13) and (14), (12) becomes 

ISl = iR'N +[-{i-f + a-f)i)' -iaV']i. (15) 

Taking the time derivative of equation (8) and the s-derivative of equation (15), and using 
equations (8) - (10), we obtain 

N' = - ^£ - = (V-'N - ijj'N) + — {fN + ^p'N) 

a a ^ 1 . (16) 

-a-f\ij\H + I- [il)V - ^V) t + I- - 1/'N) - 

2 2 - 

and 

N' = I {R'N - Ripi) - [(^7 + 07) iP" + iaV"] i - [(27 + ^7) ip' + iaV] - (^N + ^-N) . (17) 

2 

Equating (16) and (17), we obtain the following coupled equations, 

d — — 

+ a7|V'rV^ - i-^ {ipV - ipV) = -iRip - [{i-f + ^7) ip" + iaV"] (18) 

CUpi)' - ^ipip' + i^ipV' = iR' [{ij + 07) ip' + iaV] ip. (19) 

The almost impossibility of decoupling equations (18) and (19), as they are, highlights 
the enormous difficulty in fully determining the vortex motion in an analytic way from the 
Hall-Vinen equation (1). In order to make further progress, appropriate approximations are 
necessary. Recognizing that a is very small, in a first approximation, equations (18) and 
(19) may be decoupled to give the nonlinear Schrodinger equation, 

i$ ^ -|<l>P$ + 7$" (20) 
z 2 

where, 

^ = xIj + -V. (21) 

7 

Equation (20) is the same as that for the ordinary fiuid case ([E]) with the superfiuid effects 
now represented by a renormalization of the "wave function" ip, as in (21), in the first 
approximation. 
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Noting (6) and (11), (21) may be approximated as follows 

a 



1 + — (f/2+^f/3) 
7K 
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(22) 



(22) shows that, to first approximation, 

• the normal fluid velocity component along n serves to modify the curvature k; 

• the normal fluid velocity component along h serves to modify the torsion r; 
as to be expected. 

In recognition of the enormous difficulty involved in making further progress in dealing 
with equations (18) and (19) in an analytic way, we now consider a reformulation of the 
Hall-Vinen equation (1) in the extrinsic vortex filament coordinate space. 

3. Reformulation of Hall-Vinen Equation in the Extrinsic Vortex 
Filament Coordinate Space 



The extrinsic vortex filament coordinate space formulation considers only small-amplitude 
vortex motions in the LIA model. So it avoids the problems besetting the large-amplitude 
solutions due to violation of the basic premise of LIA that distant parts of the vortex filament 
remain sufficiently separated during the motion. 

Consider the vortex filament essentially aligned along the x-axis (Dmitreyev [22], Shiva- 
moggi and Heijst [23]) and take U = Uii^ + f/2i?/ + U^i./, equation (1) then becomes 

V = (1 — a') '-fut X n + at X U + a^nn — a'Uit. 

We assume the deviations from the x-axis to be small, 



r = xix + y{x, t)iy + z{x, t)iz 



We then have 



dr 

v = — = yti 



dt 

- _ c/r dr dx ^ 

ds dx ds V ^ 
We are assuming small-amplitude vortex motions, so we have 



dx 
ds 



Next, 







dx 






ds 




di 

ds 


di dx 

dx ds 


+ 


Vxx ' 


H 


+ 


z \ 





i^ + yl + zV) 



-1/2 



{jJxVxx ~l~ ZxZxx) 



1 -\{yl + zl) 



J XX ~l~ VxZxZxx 



l^-\{yl + Z^S^- [ZxVxVxx + zlZxx) 



(23) 
(24) 



(25) 
(26) 

(27) 



(28) 
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Substituting (26) - (28), equation (23) gives 

37 

2/t = - (1 - a') 7^xx + Y (l/x + ^1) ^xx + aUiZ^ - alls + a'jy^^ - a'Uiy^ (29) 

37 

zt = - {I- a') 'yy^^ + y (^^ + + <^Uiy^ + aU2 + a'jz^^ - a'Uiy^. (30) 

Putting, 

^~y + iz, V = U2 + iU3 (31) 
equations (29) and (30) can be combined to give the nonlinear Schrodinger equation 

i ($, - laV + a'U,^,) + aUi^, = (1 - a') 7$xx - ta^^.. - y (32) 

which describes the propagation of nonhnear Kelvin waves on a vortex filament in a super- 
fluid. 

Putting, 

$(x,t) = x(x,t)e^"^* (33) 

equation (32) becomes 

1 37 

- iXt + ol'UiXx) + otUiXx = (1 - a') IXxx - iot-fXxx - ylXxl^Xxx- (34) 

It may be noted that the friction term associated with a' on the left hand side in equation 
(34) can be transformed away via the Galilean transformation 

g(a;, t) ^ g(x, r), r = t — x/a'Ui (35) 

while the friction term associated with a' on the right hand side can again be eliminated 
by renormalizing the vortex strength 7. This provides a proper rationale for neglecting the 
friction term associated with a', as in [Hj, [T5] . 



4. Nonlinear Localized Structures on a Vortex Filament 

Look for a nonlinear localized Kelvin stationary wave solution 

X{x, t) = uip{x- u-fT) e^(-^-^T-)+/^- (36) 
and assume ip is slowly- varying; we then obtain from equation (34), 

(1 - ia) ^" + i {2a (1 - ia) - u]^Ij' + [c - {I - ia) - iaaUi + ifi] ^ + = 0. (37) 

We now drop ia from the factors (1 ~ia) (because a is very small compared with 1), and 
put 

a = -, (3 = a'-c, fi = — (Ui--y (38) 

The first relation implies that the velocity of propagation of the structure is twice the torsion, 
like the case with Hasimoto's [H] solution in ordinary fluid. Further, /3 is a measure of the 
curvature. Equation (37) then becomes 

r - + = (39) 
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which admits a sohtary-wave solution 



^ = W --^ sech^/|3 {x - 2(r-fT) (40) 

describing a propagating damped Kelvin kink on a vortex filament in superfluid ^He. Observe 
that the damping parameter /i vanishes when 

(i) either a = 0, ordinary fluid case, (41a) 

(ii) or u = 2Ui, special superfluid case. (41b) 

Further, the damping is symmetric with respect to the direction of propagation of the vortex 
kink (as it should be) - this symmetry is however broken by the normal fluid velocity com- 
ponent Ui along the vortex. The vortex kink growth associated with the latter aspect has 
interesting qualitative similarities with the Donnelly-Glaberson instability (Cheng et al. [24j , 
Glaberson et al. [25]) of the Kelvin waves on a vortex driven by the normal fluid velocity 
component along the vortex. 

In the second case, namely (41b), the vortex kink is undamped - the nonlinearity in the 
system, under condition (41b), balances both the dispersion and the mutual friction and 
the vortex kink amplitude and its propagation speed are determined by the normal fluid 
velocity Ui. Thus, even though the friction term associated with a, for very small a, has 
little capacity to make significant contribution to the vortex motion in a quantitative way 
it appears to be able to change vortex kink dynamics characteristics in a qualitative way. 
This feature reveals itself in a more striking way on considering the rotating planar vortex 
filament problem in superfluid "^He, as discussed below. 



5. Rotating Planar Vortex Filament 

Consider a vortex filament with shape y = y{x) lying in a plane which is rotated with 
angular velocity Q (this problem was considered by Hasimoto [2B] in the ordinary fluid case). 
We now have 

t=<l,|/,,0>— =i= (42) 

^ _ di _ didx _ / y^y^x y^^ \ 1 

""~ds dxds \ [l + ylf^^ {l + ylf'^^l ^TTyl' ^ 

Taking U = ?7iix, and using (42) and (43), equation (23) leads to 

(1 + 1/2)/ [l + yl)' 

Note that the vortex filament is rotating in the retrograde sense (i.e., opposite to that of the 
vortex core). Further, the friction term associated with a' on the left hand side can again 
be eliminated by renormalizing the vortex strength 7. 
Put, 

2/x = tan^ (45) 
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so 9 is the angle between the tangent to the vortex filament and the x-axis. Equation (44) 
then becomes 

dO 

7- aUi sin 9 = —Vly (46a) 

or 

7— -af/icos0- — + l]sin^ = 0. (46b) 
as^ as 

In the ordinary-fluid limit, equation (46b) constitutes Euler's elastica (Love [2Z]), i-e., the 
finite deformation of a plane elastic fllament of flexural rigidity B under the action of the 
thrust F applied at its ends, 

B— + FsinO = 0. (47) 

On the other hand, equation (46b) indicates a vortex growth via a normal-fluid flow 
driven instabilitjH and confirms that the friction term associated with a plays the dual roles 
of driving force and drag force ([E], [U]) and hence leads to both growth and decay of the 
vortex line length. Indeed, assuming 1^/7 1 ^ 1 and expanding in powers of 6, equation 
(46b) can be approximated by 

which is just the van der Pol equation. If the normal-fluid velocity is in the same direction as 
that of vorticity in the undisturbed vortex fllament, then there is decay/growth of the vortex 
line length if ^ 4 and vice versa if the normal-fluid velocity is in the direction opposite 
to that of the vorticity in the undisturbed vortex fllament. The friction term associated with 
a appears again to be able to change the vortex motion aspects in a qualitative way. 



6. Discussion 



Beginning with the work of Feynman [7], the phenomenological model of quantized vor- 
tices as classical vortex fllaments subject to an effective frictional force (simulating inter- 
actions with thermal excitations in superfluid ^He) has been the standard approach to in- 
vestigate vortex dynamics in superfluid ^He. The theoretical formulations developed in this 
paper provide insight into the fundamental importance of the friction term associated with 
the friction coefficient a as well as a proper rationale for the previous neglect of the other 
friction term associated with the friction coefficient a'. Further, though the friction term 
associated with a, for very small a, has little capacity to make significant contribution to the 
vortex motion in a quantitative way it appears to be able to change the vortex motion as- 
pects like the vortex kink characteristics in a qualitative way. This becomes possible because 
of the ability of the friction term to play the dual roles of driving force and drag force. 



^Experiments by Sato et al. |28| on vortex dynamics in liquid *He near indicated vortex growth due to an instability. 
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